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ABSTRACT 

We present SCUBA imaging and photometry of the eUiptical galaxy NGC 4374. 
The imaging observations are used to examine the spatial distribution of thermal 
emission from dust and the radio to infrared continuum spectrum. In the SCUBA 
850 /im image, the galaxy is found to be a point source, constraining the emission 
to less than 15" (1.5kpc for a distance of 20.74 Mpc). The simplest interpretation 
is that the 2000 fim to 850 /im SCUBA emission is synchroton from a compact core 
or inner-jet, and most unlikely thermal emission from cold diffuse dust. We cannot 
exclude free-free emission, but this would be very unexpected. The thermal emission 
from dust is shown in IRAS data and these along with the 450 /im SCUBA datum give 
a dust temperature of 35 K, corresponding to a dust mass of 1.2 x IO'^Mq. These are 
pilot observations in a program to look for cold, diffusely distributed dust in elliptical 
galaxies, however much deeper 450 fim imaging is required to investigate this in NGC 
4374. 
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1 INTRODUCTION 

The interstellar medium (ISM) in elliptical galaxies is not 
easily probed by optical observations, however observations 
at a range of other wavelengths have revealed unexpected 
amounts of gas in these galaxies (Roberts et al. 1991). With 
the gas we also expect some dust to be present, since red 
giant stars lose dust to the ISM. Small amounts of dust are 
seen in dust lanes and patches in ellipticals (Sparks et al. 
1985, Goudfrooij et al. 1994). The submillimetre (submm) 
camera SCUBA (Holland et al. 1999) provides an opportu- 
nity to test if there is cold (less than about 30 K), diffusely 
distributed dust present in elliptical galaxies as suggested 
from far-infrared observations. Dust temperatures (25 to 35 
K) and masses (~ 10* to few xIO'^Mq) were estimated from 
IRAS observations for a sample of ellipticals by Goudfrooij 
& de Jong (1995). They showed that dust masses derived 
from IRAS fluxes exceeded (by typically an order of mag- 
nitude) the dust masses estimated from optically identified 
dust lanes and patches in many ellipticals. 

However, the IRAS data are limited for studying dust 
since IRAS gave little information about the spatial distri- 
bution of the dust in galaxies, so cannot be used to check 
directly if the dust is diffusely distributed or not. We also 
note that many ellipticals contain haloes of hot gas at around 
lO'' K with masses of a few xlO^M© (Canizares, Fabbiano 
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& Trinchieri 1987). This X-ray emitting plasma is expected 
to destroy any dust grains through sputtering by hot gas 
particles in a short time (< 10^ years, Draine & Salpeter 
1979). So, for diffusely distributed dust to be present in a 
typical giant elliptical galaxy, the dust would have to be pro- 
tected or shielded from the plasma in some way. SCUBA, 
being an imaging device, offers an opportunity to map out 
the distribution of any cool dust. 

The giant elliptical NGC 4374 was chosen from the sam- 
ple of Goudfrooij & de Jong (1995) as a candidate in which 
to look for diffusely distributed dust because it has optical 
dust lanes in the central region (< 10"). The dust mass es- 
timated from the dust lanes is 3.5 x lO^M© and the dust 
mass estimated from the infrared (IRAS) fluxes (at 60 and 
100 jj.m) is 1.35 x IO^Mq (for a dust temperature of 35 K), a 
factor of ~ 4 greater than from the optically identified dust 
(Goudfrooij & de Jong 1995). Q Since most of this dust is not 
seen in the optical dust lanes, Goudfrooij & de Jong (1995) 
suggest that it must be diffusely distributed throughout the 
galaxy. Such dust would affect the colours in galaxies. How- 
ever, this affect could be difficult to disentangle from age 
and metallicity variations in the stellar population, which 



t Note that there was a typographical error in Goudfrooij & de 
Jong 1995: the dust mass they estimated from the IRAS data 
should have read LogMrf=5.13 (not 5.3) (Goudfrooij - private 
communication). We quote the correct dust mass here. 
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also produce colour changes (Worthey 1994). Therefore the 
presence of a few million solar masses of dust, distributed 
throughout a galaxy, could well have escaped optical detec- 
tion. On the other hand, spatially resolved submm obser- 
vations, sensitive to thermal emission from cool dust, may 
directly reveal the presence of any such distribution. 

NGC 4374 (M84) is an El galaxy in the Virgo cluster 
lying at a distance of 20.74 Mpc. It is classified as a steep 
spectrum Fanaroff- Riley I (FRI) radio galaxy with an in- 
frared excess. Two-sided jets emerge from its compact core, 
which is less than 2 arcseconds (Jenkins, Pooley & Riley 
1977). In the context of the Blandford & Rees (1974) jet 
model, the jets transport energy from the core to symmet- 
rical, edge-darkened lobes which extend to about 2 arcmin- 
utes and dominate the core emission (Laing & Bridle 1987, 
Wrobel 1991). The jets are perpendicular to the dust lanes. 
There is a low ionization region near the nucleus, extending 
for about 20 arcsecs along the direction of the optical dust 
lanes (Baum, Heckman & Breugel 1992). Bower et al. (1997) 
publish HST images showing the central structure and esti- 
mate a dust mass of IO^Mq (for D = 20.74 Mpc) from their 
(V-I) image. The HST images indicate central line emission 
from three components: a nuclear disk, an 'ionization cone' 
and outer filaments. 

NGC 4374 also has a hot X-ray halo with ~ IO^Mq of 
hot gas (Goudfrooij 1994). Its spectrum from the radio to 
infrared was previously studied by Knapp & Patten (1991), 
who did not detect any submm emission above that expected 
from the continuation of the power-law emission at radio 
wavelengths. They used a previous generation detector in 
the submm (Duncan et al. 1990). Knapp & Patten (1991) 
infer dust contents in a sample of radio galaxies (mostly el- 
lipticals) of the order of that found in luminous spiral galax- 
ies (~ 10® to 10* M0). They assumed a dust temperature of 
18 K (~ 15 to 20 K estimated range) for the galaxies in their 
sample and estimated a dust mass of 2 x IO^'Mq in NGC 
4374. Their low estimate of the dust temperature leads to a 
significantly larger dust mass estimate than that of Goud- 
frooij & de Jong (1995). 

In this paper we explore whether SCUBA jiggle map- 
ping at 850 and 450 yLxn and photometry at 2000, 1350, 850 
and 450 pm can confirm the presence of diffusely distributed, 
cold dust in the giant elliptical galaxy NGC 4374. We use 
the SCUBA observations together with radio and infrared 
data from the literature to place tighter constraints on the 
spectral components, including the dust temperature and 
the size of the infrared emitting region. 



2 OBSERVATIONS AND REDUCTIONS 

Imaging and photometric observations of NGC 4374 were 
obtained on 1997 September 3 as well as 1998 January 22, 
February 1, February 14 and 1999 March 19, with SCUBA, 
the Submillimetre Common-User Bolometer Array on the 
JCMT (see Table 1). In the imaging observations, we op- 
erated the 91 bolometers of the short wave array (SW) at 
450 pim and the 37 bolometers of the long wave array (LW) 
at 850 /xm, giving beam widths of 8.5" and 14.5" (FWHM) 
respectively. The arrays have a 2.3 arcminute field of view, 
and a dichroic beam-splitter allows both arrays to be used 
simultaneously. The observations employed a 64-point jiggle 



pattern, fully sampling both arrays. Similarly, the submm 
photometric observations were obtained by operating the 
central bolometers of the SW and LW arrays at 450/im and 
850/im simultaneously, employing a 9-point jiggle pattern in 
a 3 by 3 grid of 2". The 1350 /xm and 2000 /xm observations 
used the single photometry bolometers and also employed 
a 9-point jiggle pattern. Averaging the source signal in an 
area slightly larger than the beam is intended to achieve 
the best photometry accuracy under good-to-moderate see- 
ing, and, in case of the simultaneous observations with the 
LW and SW arrays, to also compensate for the very small 
offset between the arrays. During the observations, the tele- 
scope was nodded and the secondary chopped in a specified 
scheme in order to eliminate sky emission as is convention 
in submillimetre and infrared astronomy. 

Residual sky emission was removed off-line by using 
quiet SCUBA array bolometers in which there was no source 
emission, usually the bolometers in the first ring (for LW) 
and second ring (for SW) from the centre (Jenness et al. 
1998, Robson et al. 1998). The pointing stability was checked 
before and after each map and before each photometry ob- 
servation. The focus was checked every three hours or when 
there were noticeable dome temperature fluctuations. The 
atmospheric opacity was measured regularly by performing 
SCUBA skydips. The sky monitor at 225 GHz on the Cal- 
tech Submillimetre Observatory (commonly known as CSO 
tau or rcso, see Masson 1993) updates every 15 minutes 
and so was used to monitor sudden changes in atmospheric 
opacity, and, in cases where the SCUBA skydips were not 
measured or produced poor fits, the CSO data were used to 
extrapolate SCUBA atmospheric opacities as described in 
the SCUBA observing manual. 

The dedicated SCUBA data reduction software SURF 
(Jenness et al. 1998), as well as the Starlink packages Kappa 
and Figaro were used to reduce and analyse the observa- 
tions. The data reduction consisted of first fiat-fielding, cor- 
recting for atmospheric extinction, and then residual sky 
emission removal (see above). Noisy bolometers and inte- 
grations were blanked and spikes removed. For the 1350 /im 
and 2000 /im photometry observations no residual sky emis- 
sion removal was done as the observation used only a sin- 
gle bolometer. Spike removal was performed by clipping the 
data at a specified sigma. The resulting data were calibrated 
using instrumental gains that were determined from beam 
maps of Mars or Uranus nightly and in the same observa- 
tion mode as the target observation (see Table 1). Plane- 
tary fiuxes for each filter were obtained using the JCMT 
utility program FLUXES. During the January and Febru- 
ary 1998 runs, the planets were not available and the JCMT 
secondary calibrators CRL2688 and IRC-l-10216 were used 
(SandeU 1994; Sandell in preparation). On nights when the 
Tcso > 0.07 (see above and Table 1) the atmosphere was 
opaque to 450/im emission and observations could be made 
only at longer wavelengths. 



3 RESULTS AND DISCUSSION 

Figure 1 shows the 850 fj,ra image of NGC 4374 together with 
the beam-width determined from observations of 3C 273, a 
standard point source for the JCMT, made in the same way 
and conditions as NGC 4374. The rms uncertainty on the 
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Table 1. Submillimetre fluxes (Janskys) for NGC 4374 from SCUBA 
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Figure 1. The SCUBA 850 ^tm image of NGC 4374. The observation was made on 1998 February 1 when rcso was ~ 0.05. 



map is less than 20mJy/beam. The image of NGC 4374 
has small northeastern and southwestern spurs that are not 
real features of the galaxy but a manifestation of the beam 
smearing as a result of chopping. These spurs are in the 
chop direction and are detected at a much lower level than 
the central emission. An examination of the point spread 
function of the NGC 4374 image yielded FWHMminor-axis = 
15.0" and FWHMmajor-axis = 15.7" The result shows that 
the image of NGC 4374 is no more than 1" greater than the 
JCMT beam. This insignificant spatial extension is probably 
a result of a small pointing drift during the long observation 
(~ Ihr). We conclude that NGC 4374 is not extended in the 
SCUBA 850/im image which constrains the spatial origin of 
the observed emission to less than 15" (1.5 kpc for a distance 
of 20.74 Mpc) in diameter. It is worth noting that although 



the SCUBA beam is ~ 15" at 850 /im, the arrays have a 
2.3 arcminute field of view, and therefore, at the distance 
of NGC 4374, we do not detect any diffused emission to 
~ 13 kpc in extent, as we observed with the array pointed at 
the centre of the galaxy and our observations were spatially 
sensitive to the entire field of view of SCUBA. 

The results from the imaging and photometric observa- 
tions are shown in Table 1. The uncertainties are a quadratic 
sum of the uncertainty arising from the measured signal-to- 
noise-ratio and a systematic calibration uncertainty, which 
varies from 10% at 850 [xm and lower, to 27% at 450 /im. The 
450 p-tcl fiuxes listed for the 1998 Feb 01 and 14 and Jan 22 
runs are 3(t upper limits. The submillimetre flux variation 
over the four observing runs was within the errors on the 
fiuxes. 
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3.1 Decomposing the SED 

Figure 2 shows the spectral energy distribution (SED) for 
NGC 4374, including data from radio (Jenkins, Pooley & 
Riley 1977) to infrared (IRAS) wavelengths. The infrared 
fluxes were obtained from NED (the NASA/IPAC Extra- 
galactic Database) and have been interpreted as thermal 
emission from dust (Knapp & Patten 1991, Goudfrooij 
1994). At longer radio wavelengths, the lobes dominate the 
integrated flux from NGC 4374 (by a factor of 8 at 5 GHz, 
Wrobel 1991) and give an excellent power law of spectral 
index (S^ cx of Ogi'^GHz" = -0.6 ± 0.03, as shown by 
the dotted line in Figure 2. This spectral index is consistent 
with the classification of NGC 4374 as a steep spectrum 
radio galaxy. The good fit suggests variability larger than 
the errors on the fluxes is uncommon for the synchrotron 
power-law of the radio emission. 

The submm data points for fluxes within the SCUBA 
beam and previous submm measurements fall below an ex- 
trapolation of this steep radio power-law from the integrated 
flux. We see no sign of extended emission from the radio 
lobes on any of our maps, and the 3a upper limit on the 
850 ^m integrated flux within a region approximately the 
size of the lobes is 0.6 Jy. This is much higher than the 
extrapolation of the radio power-law from the integrated 
fluxes, showing that we are insensitive to the radio lobes 
(Figure 1). As we show later, this indicates that if the ob- 
served SCUBA 850 fj,m emission is non-thermal in nature it 
is dominated by the compact core (or inner-jet). 



3.1.1 The Core SED 

In Figure 2 we have plotted the single available core radio 
flux for NGC 4374, at 5GHz (Jenkins, Pooley & Riley 1977, 
Wrobel 1991); the circles are the means of the SCUBA data 
listed in Table 1. The 5 GHz core flux together with the 
SCUBA data give a power-law slope of a = —0.23 ± 0.03. 
In compact radio galaxies, such flat spectra have been found 
to comprise spectra from many unresolved components, each 
with its own synchrotron spectrum, the sum of which results 
in the observed flat power-law. Although we do not wish to 
over-interpret our results, since we only have one radio core 
point and the SCUBA data, we note that the observed flat- 
ter power-law of the NGC 4374 core could be a consequence 
of such unresolved synchrotron components. At 5GHz, the 
northern jet in NGC 4374 is brighter than the southern jet 
for the flrst 10 arcseconds (Bridle and Perley 1984). This 
one-sided inner-jet may comprise such unresolved compo- 
nents. 

It is worth noting that the SCUBA fluxes at 850 fim, 
1350 pm and 2000 /im alone give a spectral index of 
QiteGHz ~ —0.16 ± 0.32 in the submm, forming a flat spec- 
trum resembling that of a free-free emission component. 
While this 2000 /xm to 850 /xm spectral index is consistent 
with that of free-free emission (a = —0.1), the large uncer- 
tainty on its value limits our ability to conclusively prove 
its origin. Our 450 /im point lies well below this a = —0.16 
spectrum, which does not help to prove or disprove its origin, 
but instead adds the complication that this free-free com- 
ponent might have a high-frequency cut-off between 850 
and 450 pm. In steep spectrum radio galaxies such as NGC 
4374, it is uncharacteristic for a free-free component to dom- 



Table 2. Parameters determined from the model fit to the SED 
for NGC 4374 

Parameter Value Error Units 



a -0.23 ±0.03 - 

n 3.96E-11 - Steradians 

1.46 - arcseconds 

T 35 ±5 K 

inate the power-law in the mm-submm wavelength range. 
Therefore if this component is real it would have to come 
from very high ionization and it would be an interesting dis- 
covery. We cross-correlated the 'submm flat spectrum com- 
ponent' in Figure 2 with the well known spectrum of 3C273 
using data taken with SCUBA on 1998 February 15, just 
one night after and under very similar conditions as the 
photometry observations of NGC 4374. The spectrum of 
3C273 gives a much steeper slope and a spectral index of 
Qi46GHz ~ — 0.7 ± 0.1, as expected of this variable object 
during its quiescence phase (Robson et al. 1993). Therefore 
with our current, limited data we cannot say the 850 /im, 
1350 fim and 2000 /im fluxes come from significant free-free 
emission, if any at all. 

To model the radio to infrared SED we used a combi- 
nation of a power-law (non-thermal radio emission from the 
active nucleus) plus greybody (re-processed emission from 
dust in the infrared, Hughes, Gear & Robson 1994). The 
greybody allows for the fact that we see to different physi- 
cal depths at different wavelengths, for a given optical depth. 
The composite model fiux is then: 



(1) 



In equation (1) C is the normalisation for the power- law 
component, a the power-law spectral slope, fl the solid an- 
gle for the greybody component, Bv{T) the Planck function 
at temperature T, Ao the wavelength at which the optical 
depth is unity and /3 the emissivity index of the grains. The 
normalisation factors (C and Q.) are obtained by forcing the 
model to agree with the 5GHz core radio and IR (60 and 
100 jUm) fluxes respectively, since the two components dom- 
inate in these different wavebands. The power-law index and 
its uncertainty were estimated from linear regression fits to 
the core radio plus SCUBA data in a log-log plot. We assume 
that the IRAS emission comes from a region less than ~ 10" 
due to lack of extension seen at 450 /xm. The best fit tem- 
perature and its uncertainty for the greybody component 
were estimated by eye from plots of the composite model 
versus the data. Other parameter values (Ao = 7.9 pm and 
P = 1.3 ~ ±0.5) were fixed at the values given in Hughes 
et al. (1994) in their study of M82. The parameter values 
determined from the model are summarised in Table 2, and 
the model is plotted as a solid line in Figure 2. 

We note that the greybody fit demands that the warm 
dust emission comes from a very compact zone, much smaller 
that the SCUBA 450 /xm resolution (see Table 2). Because 
the 450 yLxn data point lies below the sum of the thermal 
and non-thermal emission (see Figure 2) , the simplest expla- 
nation is that the synchrotron spectrum steepens between 
850 pm and 450 /xm (or somewhat long-ward of 850 /im, 
given the uncertainties) . Also, some small part of the 450 /xm 
emission could be due to low surface brightness dust of 
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Figure 2. Spectral energy distribution from radio to infrared for the core of the elliptical galaxy NGC 4374. 
Filled circles are our data from SCUBA observations and the sources of other data are described in the text. 
The solid line is a composite model of a power-law {a = —0.23) plus (35 K) greybody. This model fits the 
data quite well for the radio through to 60/im infrared detections. The dashed line shows the 35K greybody 
alone, showing that the SCUBA data lie above this component. The dotted line is a power law of slope —0.6, 
passing through the integrated radio fluxes to 31 GHz. The downward arrow indicates a 3 sigma upper limit 
at 25fim. 



~ 15K to 20 K. Better spectral coverage and higher signal- 
to-noise data in the submm are clearly needed to test more 
complex spectral models. 



3.2 Dust Mass and its Implications 

The temperature that fits the IRAS data and is constrained 
by the SCUBA 450 /im measurement is ~ 35 ± 5 K (see Fig- 
ure 2). Not surprisingly, this is the same temperature as 
Goudfrooij (1994) found for this galaxy from his analysis 
of the IRAS data alone. Changing /3 between 2.0 and 1.0 
makes very little difference to the best fit temperature. The 
mass of emitting dust Md can be derived from a simple 
model adapted from Hildebrand (1983), where 



S,D' 



(2) 



In equation (2) Si, is the measured fiux, D the distance to the 
source (20.74 Mpc for NGC 4374 as assumed by Goudfrooij 
1994), B^iT) the Planck function and fed the grain mass 
absorption coefficient. We assumed k^*^ ^'^ — 2.5m^kg~^ 
(Hildebrand 1983) and estimated a dust mass of 1.2xlO^M0 
for Sioo Mm = 0.98±0.21 Jy and T = 35 K. The dust mass we 
have calculated for NGC 4374 is about the same as Goud- 
frooij & de Jong (1995), as expected since they used the 
same temperature. Furthermore, it is similar to that found 
in a recent HST {V - 1) study of dust lane in NGC 4374 by 
Bower et al. (1997), who found the dust mass to be in agree- 
ment with that derived from IRAS data, in contrast with 
the dust deficit found in other optical studies of ellipticals 
(Goudfrooij & de Jong 1995). However for NGC 4374 it is 
at least one order of magnitude lower than Knapp & Patten 
(1991), who assumed a cooler temperature of 18 K (based 
on objects in their sample of nearby radio galaxies for which 
1350 /im, 800 /xm and 450 /xm fluxes were detected). Also, 
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they obtained the higher mass even though they assumed a 
smaller distance of 13 Mpc. While the low dust temperatures 
estimated by Knapp & Patten (15 to 20 K) are similar to 
the Galactic Plane value of T = 19 K, they are inconsistent 
with the SCUBA plus IRAS observations of NGC 4374 if 
the dust is assumed to be all a single temperature. 

As for the extended radio lobe emission, we have not 
detected the extended low-level emission from diffusely 
distributed dust with SCUBA. We have shown that the 
2000 to 850 /im data are most unlikely due to diffuse cold 
dust, but some small part of the 450 ^.m emission could be 
due to dust of ~ 15 K to 20 K. Also, our 3(j upper limit to the 
surface brightness at 450 /im is too high and not useful for 
setting dust mass upper limits on this possibly colder dust. 
Deeper SCUBA imaging observations are clearly needed to 
address this question and will be undertaken in the 2000 
observing season. 



4 CONCLUSIONS 

Following the suggestion that elliptical galaxies may contain 
diffusely distributed dust (Goudfrooij 1994) we searched for 
this dust with submm imaging observations of the elliptical 
galaxy NGC 4374, using SCUBA on the JCMT. We have not 
detected low-level, diffusely distributed dust with SCUBA. 
The emission at 850 /im is spatially unresolved (diameter 
< 15"; 1.5kpc). 

Adding the SCUBA submm data to existing radio 
through IR data for this galaxy we can constrain the dust 
component to a single temperature of ~ 30 to 40 K, imply- 
ing a dust mass of ~ 1.2 x lO^M©. The model fitting in 
section 3.1 gives an angular extent of the dust to be ~ 1.5" 
in diameter. Mindful that this result is too simplistic for a 
radiative model of an AGN torus, we note in passing that if 
we assume a co-mixing of the molecular gas and dust, this 
constrains the size of any molecular torus around the AGN 
core of NGC 4374 to ~ 150 pc in diameter. 

The mm-submm observations show a flat spectral index 
that is consistent with the 5 GHz radio core flux. The spec- 
tral index of the 850 jim to 2000 /im fluxes alone is consistent 
with that expected from free-free emission. The possibility of 
a free-free component is very unusual and intriguing. Given 
the care we took over calibration, we see no systematic rea- 
sons for our peculiar mm-submm SED. Future observations 
of NGC 4374 will aim to achieve very high signal-to-noise 
observations in order to determine the spectral index to a 
much higher accuracy and therefore prove or disprove the 
presence of a third spectral component. 

In future observations of dusty ellipticals, we will con- 
centrate on those with luminous infrared emission, coldest 
dust and less non-thermal radio emission in order to try to 
resolve the distribution of cold dust in ellipticals. We will 
obtain short mm images of the central regions of NGC 4374 
to determine the spatial extent of the emitting region(s) and 
deeper SCUBA imaging observations to address the possi- 
bility of extended low- level emission from cold dust. These 
observations together with ISO data will be used in or- 
der to get a tighter constraint, and thus a better handle, 
of the radio to submm SED. These observations have im- 
portant consequences for the interpretation of colours and 
colour gradients in elliptical galaxies, which are attributed 



to age and metallicity changes in the absence of diffusely 
distributed dust. 
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